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ABSTRACT 

Populations of five related species belonging to the Gam - 
marus pulex group from the Iberian Peninsula and southern 
France have been studied electrophoretically at 21 enzyme 
loci. Morphologically distinct forms from the same side of 
the Pyrenees proved to be genetically more similar than 
morphologically very similar forms from opposite sides of 
this barrier. 

With reference to published data, genetic differentiation 
between French and Iberian species was observed at the 
generic level. Inter-areal intraspecific genetic differentia¬ 
tion in Iberian species may warrant recognition of 
subspecies or sibling species. However, cross-breeding 
experiments performed did not assess levels of reproduc¬ 
tive isolation of the forms involved. 

A dendrogram derived from genetic distances ( D ) was 
compared with geological and paleoclimatological 
evidence. Although hypotheses derived from the alter¬ 
native data bases do not conflict, they give little informa¬ 
tion with regard to evolutionary rates of the loci involved. 

RESUME 

Des populations de cinq especes apparentees appartenant 
au groupe Gammarus pulex provenant de la peninsule Iberi- 
que et du Midi de la France ont ete etudiees par elec- 
trophorese a 21 loci d’enzymes. Des formes morphologi- 
quement distinctes du meme versant des Pyrenees se sont 
montrees genetiquement plus proches que des formes 
morphologiquement fort semblantes du versant oppose de 
cette barriere. 

En se basant sur des donnees publiees, la differentiation 
genetique entre les especes iberiques et fran^aises a ete 
constatee au niveau generique. La differentiation geneti¬ 
que intra-specifique entre aires de repartition dans le cas 
des especes iberiques, peut permettre de reconnaitre des 
sous-especes ou des especes-jumelles. Cependant, des 
essais d’hybridation effectuees n’ont pu etablir une isola¬ 
tion reproductive, a divers niveaux, des formes con- 
cernees. 

Un dendrogramme deriv6 de distances genetiques (D) a 
ete confronte avec des donnees geologiques et paleoclima- 


tologiques. Bien que les hypotheses basees sur d’autres 
series de donnees diverses ne se contredisent pas, elles 
apportent peu d’information sur les rapports evolutifs des 
loci concernes. 


INTRODUCTION 

The Gammarus pulex group (sensu Karaman & 
Pinkster, 1977) of the Iberian Peninsula is rep¬ 
resented by G. gauthieri S. Karaman, 1935 and 
G. ibericus Margalef, 1951. G. gauthieri is known 
from Tunisia, Algeria, Morocco and many 
more or less isolated localities in Spain 
(Pinkster, 1971). Until 1974 G. ibericus was 
known only from the type locality. In that year, 
the species was recorded from many localities in 
the Massif Central, France (Goedmakers, 
1974). Indeed, individuals from these popula¬ 
tions resembled G. ibericus in almost every mor¬ 
phological detail. Pinkster & Scholl (1984), how¬ 
ever, demonstrated that a sample of individuals 
from one of these localities (station 2 in the 
present study; table I, fig. 1A) was genetically 
distinct and reproductively isolated from G. 
ibericus . This form was described as G. orinos 
Pinkster & Scholl, 1984. 

Geologically, the Iberian Peninsula and the 
rest of the European continent have developed 
quite independently (e.g. Vandenberg, 1979; 
Dercourt et al., 1985). Moreover, it is well 
known that there exists a rather close relation¬ 
ship between the freshwater fauna of the 
Iberian Peninsula and that of North Africa 
(Margalef, 1983; Pinkster, 1971; Garcia de 
Jalon Lastra & Gonzalez del Tanago, 1986). 
Reproductively isolated populations of the G. 
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pulex group in France (Scheepmaker, un¬ 
published) differ genetically much less from 
each other than G. orinos and G. ibericus 
(Pinkster & Scholl, 1984). The possibility exists 
that morphologically rather different forms 
from the Iberian Peninsula (e.g. G. ibericus and 
G. gauthieri) may be genetically more similar 
than even morphologically more comparable 
forms on either (the Iberian and the European) 
side of the Pyrenees (i.e. G. ibericus versus G. 
orinos , G. stupendus , and a morphologically 
distinct form restricted to the Pyrenees which is 
referred to as G. cf. Jossarum by Goedmakers & 
Roux, 1975). 

The primary aim of this study was to deter¬ 
mine if the Iberian members of the G. pulex 
group form a genetic entity separate from those 
of the rest of the European continent. Secondly, 
as some members of the G. pulex group, 
although genetically homogeneous, show con¬ 
siderable morphological variation and vice 
versa (Scheepmaker, 1987), the possibility of 
geographically proximate populations of G. 
ibericus and G. gauthieri being closer related to 
each other than allopatric populations of G. 
gauthieri was investigated. 

Genetic variation at 21 presumptive gene loci 
was investigated using starch gel enzyme elec¬ 
trophoresis. In order to evaluate the biological 
significance of the genetic variation recorded, 
cross-breeding experiments have been carried 
out between two geographically distanct pop¬ 
ulations of G. gauthieri , two genetically rather 
different populations of G. ibericus and two more 
or less sympatric populations of G. gauthieri and 
G. ibericus. From these data the relationships 
among the samples studied were evaluated. 

MATERIAL AND METHODS 
Sampling and collection sites 

Sampling for electrophoretic studies and cross-breeding 
experiments was carried out based on the procedures of 
Scheepmaker (1987). The collection sites and sampling 
dates are listed in table I and fig. 1. Sampling localities for 
electrophoretic analysis have been chosen according to 
Pinkster (1971), Goedmakers (1974) and Goedmakers & 
Roux (1975), who did extensive sampling on the Iberian 
Peninsula, in the Massif Central (central France) and in 


the northern Pyrenees (southern France). In the follow¬ 
ing, samples of individuals of populations will be referred 
to as “populations”. 

Morphologically similar forms of G. gauthieri were 
sampled in the Valladolid area (fig. 1C: stations 10, 11 
and 12) and the Serranfa de Cuenca (fig. ID: stations 6 
and 7). Morphologically, populations from these areas dif¬ 
fer slightly in the setation of P5-P7 (figs. 2D, E). Except 
for the length of the inner ramus of uropod 3 (fig. 2N), 
populations from both areas resemble G. jossarum Koch in 
many aspects. 

G. ibericus differs from G. gauthieri by the setosity of P5- 
P7 (figs. 2D, E, F) and the absence of setae on the dorsal 
side of the telson (figs. 2Q, S). Morphologically identical 
populations of G. ibericus were sampled in different 
drainage systems in the Serranfa de Cuenca (table I, fig. 
ID: stations 8 and 9). G. orinos (table I, fig. 1 A: station 2), 
G. cf. jossarum (table I, fig. IB: stations 3, 4, 5) and G. 
stupendus (table I, fig. 1A: station 1) are morphologically 
related forms of which above all G. orinos resembles G. 
ibericus in many details (e.g. figs. 2G, F, M, O). Com¬ 
parison of the differentiation among these forms as for 
morphological characters is summarized in table II. 

Electrophoresis 

The genetic variation at 21 enzyme loci was studied using 
both horizontal and vertical starch gel electrophoresis. 
Electrophoresis generally followed Siegismund et al., 1985 
(horizontal starch gel electrophoresis) and Bulnheim & 
Scholl, 1981 (vertical starch gel electrophoresis). Starch 
gels (Connaught 12% w/v) were run for 15-17 hours at 20 
V cm’ 1 (buffer IV; 8 V cm' 1 for vertical electrophoresis), 
5 to 6 V cm' 1 (buffers I, II, III, V and VI) and 3 to 4 V 
cm' 1 (buffer VII). The enzymes assayed are enumerated 
in table III along with the corresponding buffer systems. 
Staining techniques of ADA, ALP, APK, EST, GOT, 
GPI, LAP, MDH, ME, MPI, and PEP in general fol¬ 
lowed Siegismund et al. (1985). The segregation of the 
Est-1 and Est-2 loci requires long pathways. Therefore, in 
horizontal electrophoresis, extra long molds were used. 
Samples from station 1 serve as “marker” in all our 
studies. Thus, sample sizes of station 1 are nearly always 
exceeding 99 (table IV). 

Individuals in some samples gave very poor results for 
the GDH, MDH and ME stainings, thus the substrate 
concentration was generally doubled. For the PEP stain¬ 
ings, 1-leucyl-glycyl-glycine and 1-leucyl-glycyl-tyrosine 
were used as substrate. GDH, HK, 6PGD and PK stain¬ 
ings were according to Shaw & Prasad (1970), Jelnes 
(1971), Giblett (1969) and modified after Brewer (1970), 
respectively. 

Analysis oj allozyme variation 

Genetic interpretation was inferential. Electromorph fre¬ 
quencies and matrices of D and / (according to Nei, 1972) 
and of Rogers’ (1972) genetic distance were calculated 
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Table I 


Sampling localities and species distribution. 


Station 

no. 

Species 

Country 

Prov. 

/dept. 

Drainage 

system 

Locality 

description 

Sampling 

date 

1 

G. stupendus 

France 

Var 

Gapeau 

Source du Gapeau, alongside road 
D2, 2 km S. of Signes, 40 km N. of 
Toulon 

29-VII-’86 

2 

G. orinos 

France 

Ardeche 

Ardeche 

Rivulet draining into the Ardeche, at 
crossing with road D533, 4 km N.W. 
of Vals-les-Bains 

10-VII-’86 

7-V-’86 

3 

G. cf. Jossarum 

France 

Ariege 

Ariege 

Confluent of the Ariege, 0.5 km N. of 
Alliat, ca. 10 km W. of Tarascon 

7-VIII-’87 

16-X-’87 

4 

G. cf. Jossarum 

France 

Ariege 

Ariege 

Tricklet alongside road D18, 0.5 km 
W. of Sue, ca. 20 km W. of Tarascon 

7-V-’86 

7-VIII-’87 

5 

G. cf. Jossarum 

France 

Ariege 

Garonne 

Rivulet crossing road D18, 1.2 km E. 
of Aulus-les-Bains, ca. 30 km W. of 
Tarascon. 

7-V-’86 

6 

G. gauthieri 

Spain 

Cuenca 

Tajo 

Rio Albalate, alongside road N320 at 
Arralba, ca. 32 km N.N.W. of 
Cuenca. 

9-V-’86 

7 

G. gauthieri 

Spain 

Cuenca 

Guadiana 

Rio Cigiiela, alongside road N400 at 
Naharros, ca. 35 km W. of Cuenca. 

10-V-’86 

29-V-’86 

8 

G. ibericus 

Spain 

Cuenca 

Jucar 

Rio Jucar, at its confluence with the 
Laguna de Una, ca. 36 km N. of 
Cuenca 

10-V-’86 

29-IV-’87 

9 

G. ibericus 

Spain 

Cuenca 

Tajo 

Spring called Canada de las Tablas, 
draining into the Rio Cuerbo, 12 km 
N. ofTragacete 

28-IV-’87 

10 

G. gauthieri 

Spain 

Valladolid 

Duero 

Rio Cega, crossing road Cl 12 near 
Iscar, ca. 12 km W. of Cuellar 

13-V-’86 

29-IV-’87 

11 

G. gauthieri 

Spain 

Segovia 

Duero 

Rio Cerquilla, crossing road C601, 2 
km S. of Cuellar 

12-V-’86 

12 

G. gauthieri 

Spain 

Valladolid 

Duero 

Confluent of Rio Pisuerga, crossing 
road Va901, 2 km S. of Cigales 

13-V-’86 


Table II 


Morphological differentiation of some taxonomically important characters in G. orinos, G. cf. Jossarum , G. stupendus , 

G. ibericus, and G. gauthieri. 


Station 

G. orinos 

2 

G. cf. Jossarum 
3. 4, 5 

G. stupendus 

1 

G. ibericus 

8, 9 

G. gauthieri 
6, 7 

G. gauthieri 

10, 11, 12 

Setosity A2 

moderate 

moderate 

moderate 

moderate 

little 

little 

Calceoli 

+ 

— 

— 

— 

+ 

+ 

Setosity P5-7 

long 

moderate 

long 

long 

absent 

little 

Setae longer than spines 

+ + 

+ 

+ + 

+ + 

— 

— 

Uropod 3: endop./exop. 

1/2 

1/2 

1/2 

2/3-3Z4 

3/4 

3/4 

Setae on dorsal surface telson 

+ 

+ 

+ 

— 

+ 

+ 
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Fig. 1. A-D, Study areas and sampling localities (cf. table I). 


with the computer program BIOSYS-1 (Swofford & 
Selander, 1981). From these data, an UPGMA den¬ 
drogram (Sneath & Sokal, 1973) of Nei’s genetic distance 
(D) (Nei, 1972) and a distance Wagner network (Farris, 
1972) based upon Rogers’ genetic distance (Rogers, 1972) 
were constructed. A “character state” transformation 
scheme was calculated with the aid of the program Jelly 
(Ellis, 1987). A chi-square test was performed to test 
observed genotype frequencies and those expected under 
Hardy-Weinberg equilibrium. Electromorph frequencies 
of samples taken at stations 3, 4, 7, 8 and 10 (table I, figs. 
1B-D) did not differ significantly in 1986 and 1987, thus 
these data are combined in the analysis. 

Nei’s D has often been criticized because of its non- 
metricity which makes it unsuitable for deriving phylo¬ 
genetic trees. Moreover, both Nei’s distance coefficient 
and the UPGMA method require the assumption of 
homogeneous evolutionary rates (see Farris, 1981; Swof¬ 
ford, 1981; Mickevich & Miller, 1981; Buth, 1984; 
Rogers, 1986). Nevertheless, Nei’s D- and /-values 


(whether or not in combination with the UPGMA 
method) have been widely used by authors for studying 
speciation and for estimating divergence time (e.g. 
Thorpe, 1982; for Crustacea, among others, Sbordoni et 
al., 1980; Siegismund et al., 1985; Bert, 1986). 

The Wagner distance procedure maximizes the infor¬ 
mation contained in the matrix of genetic distances and 
does not depend on the assumption of homogeneous 
evolutionary rates. This method results in the formation 
of the most parsimonious trees, i.e., trees with a minimum 
length. 

Shortcomings of methods based upon distance matrices 
are the generation of “impossible ancestors” (Farris, 
1981; Mickevich & Mitter, 1981; Swofford & Berlocher, 
1987) and the neglection of less common electromorphs: a 
group of taxa can exhibit uniformly high pairwise genetic 
similarities and still possess a very different combination 
of electromorphs (Swofford & Berlocher, 1987). 

The Jelly program (Ellis, 1987) copes with these 
problems. This method, comparable to the ROGAD 
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Fig. 2. Morphological differentiation of some taxonomically important characters in males of G. orinos (station 2), G. 
cf . fossarum (station 3), G. stupendus (station 1), G. ibericus (station 8) and G. gauthieri (stations 6+ 11). 

A-F: pereiopod 7 of G. stupendus , G. cf. fossarum, G. orinos , G. gauthieri from station 11, G. gauthieri from station 6, and 
G. ibericus , respectively; 

G-K: antenna 2 of G. cf. fossarum , G. ibericus , G. orinos, G. gauthieri from station 6, and G. stupendus , respectively; 
L-P: uropod 3 of G. cf. fossarum , G. ibericus , G. gauthieri from station 6, G. orinos and G. stupendus , respectively; 
Q-U: telson of G. gauthieri from station 6, G. cf. fossarum, G. ibericus , G. stupendus and G. orinos , respectively. 
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Table III 


Enzymes studied and buffer systems employed. 


Enzyme 

Abbreviation 

E.C. no. 

Buffer system* 

Number of loci 

Adenosine deaminase 

ADA 

3.5.4.4 

IV 

2 

Alkaline phosphatase 

ALP 

3.1.3.1 

VI 

2 

Arginine phosphate kinase 

APK 

2.7.3.3 

II, IV, V 

1 

Esterase 

EST 

3.1.1.1 

III 

2 

Glutamate dehydrogenase 

GDH 

1.1.1.47 

III 

1 

Glutamic oxaloacetic transaminase 

GOT 

2.6.1.1 

I 

2 

Glucose phosphate isomerase 

GPI 

5.3.1.11 

II, V 

1 

Hexokinase 

HK 

2.7.1.1 

IV 

3 

Leucine aminopeptidase 

LAP 

3.4.1.1 

IV 

1 

Malate dehydrogenase 

MDH 

1.1.1.37 

I 

2 

Malic enzyme 

ME 

1.1.1.40 

II, IV 

1 

Mannose phosphate isomerase 

MPI 

5.3.1.8 

II 

1 

Peptidase 

PEP 

3.4.11/13 

V 

5 

6 Phosphogluconate dehydrogenase 

6PGD 

1.1.1.44 

VII 

1 

Pyruvate kinase 

PK 

2.7.1.40 

II, V, VII 

1 


* The following buffer systems were used: (I) Amine citrate buffer pH 6.1 (Clayton & Tretiak, 1972); (II) Tris maleate 
buffer pH 7.2 (Harris & Hopkinson, 1976); (III) Tris citrate buffer pH 8.0 (Selander et al., 1971); (IV) Tris EDTA 
borate buffer pH 8.9 (Ayala et al., 1972; modification after Scholl et al., 1978); (V) phosphate buffer pH 6.7 
(Selander et al., 1971); (VI) Tris citrate borate discontinuous buffer system pH 8.2-8.7 (Poulik, 1957); (VII) 
Histidine sodium citrate buffer (Brewer, 1970), with the addition of NADP (5 mg/1 in gel, 10 mg/1 in cathodal tray 
buffer). 


method (Rogers, 1984; Swofford & Berlocher, 1987), 
generates Wagner character states from elecromorph fre¬ 
quency characters using the HAP algorithm (Rogers, 
1984). It considers frequencies of a particular elec- 
tromorph as a character and any combination of elec- 
tromorph frequencies as a particular character state. The 
frequency distributions (“character states”) attributed to 
each of the HTU’s (Hypothetical Taxonomic Units) are 
not selected from one of the “states” occurring in the 
OTU’s (Operational Taxonomic Units), but have unique 
values of their own. These “character states” consist of 
arrays of electromorph frequencies in which the frequen¬ 
cies sum up to one, thus satisfying the “additivity require¬ 
ment” (see Swofford & Berlocher, 1987). This method has 
also the advantages of the Hennigian method as proposed 
by Patton & Avise (1983), but it copes with criticisms by 
Swofford & Berlocher (1987) by considering all characters 
simultaneously. 

The Jelly program starts by generating a user-defined 
number of networks, which are stocked in a so-called 
“arena”. The concept of the arena is based on the model 
of Darwinian selection. The arena is an imaginary zone of 
memory, where a population of competing networks 
resides. Then, in a user-defined number of cycles, one of 
the networks is selected; the chance of being selected is 
maximal for the best network. Part of the OTU’s is 
removed from the network and then refitted in an attempt 


to find a better network. The resulting network replaces 
the worst network in the arena. When no further improve¬ 
ment (viz. decrease in length) occurs in the residing net¬ 
works, the program is stopped. The results obtained by 
the three methods discussed above are compared in order 
to evaluate differences due to the application of a par¬ 
ticular method. 

Cross-breeding experiments 

Cross-breeding experiments followed Pinkster (1983) and 
were carried out (a) between morphologically similar but 
geographically distant G. gauthieri populations from the 
Valladolid area (station 10; fig. 1C) and the Serrama de 
Cuenca (station 7; fig. ID); (b) between morphologically 
different but geographically neighbouring populations of 
G. ibericus (station 8; fig. ID) and G. gauthieri (station 7; 
fig. ID) and (c) between morphologically identical but 
genetically different populations of G . ibericus (stations 8 
and 9; fig. ID). 


RESULTS 

Twenty-one gene loci from seventeen enzyme 
systems were scorable in all populations. Elec- 
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tromorph frequencies and relative mobilities of 
the enzymes coded for by these loci are given in 
table IV. Electromorphs of the Lap locus were 
monomorphic for all populations. Two or more 
loci were resolved in eight of the enzyme 
systems: ADA, ALP, EST, GOT, HK, MDH 


and both PEP stainings (table III). Among 
these, the locus that codes for the enzyme with 
the highest anodal mobility has been numbered 
1. The gene products for two enzyme systems 
moved cathodically: GOT-2 and MDH-2. 
Enzyme products of the Pep loci have been 


Table IV 

Electromorph frequency distribution at 20 enzyme loci. N = sample size; h - heterozygosity per locus (direct count); 
H = mean heterozygosity over all (including monomorphic) loci; a = fastest moving electromorph; * = signifant depar¬ 
ture from Hardy-Weinberg distribution (P<0.05). 


, Station no. 

Locus - 


electromorph 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Ada-1 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

1.00 

1.00 

1.00 


N 

99 

34 

20 

37 

27 

25 

28 

47 

32 

30 

47 

34 

Alp-1 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 


N 

99 

47 

50 

50 

40 

45 

33 

28 

26 

45 

59 

48 

Alp-2 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

0.61 

0.70 

1.00 

1.00 

1.00 

1.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.39 

0.30 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.39 

0.42 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

99 

47 

50 

50 

40 

41 

33 

59 

45 

45 

59 

48 

Apk 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.96 

0.96 

0.97 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.04 

0.04 

0.03 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.08 

0.08 

0.05 


N 

99 

50 

45 

53 

42 

57 

38 

81 

33 

39 

51 

39 

Est-1 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 


c 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

99 

8 

16 

27 

33 

18 

28 

51 

51 

38 

22 

17 

Est-2 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.00 

1.00 

1.00 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 


N 

99 

8 

16 

27 

33 

18 

28 

51 

51 

38 

22 

17 

Gdh 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 


N 

99 

71 

58 

59 

94 

82 

57 

80 

57 

70 

63 

33 

Got-1 

a 

0.00 

0.01 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.00 

1.00 

1.00 

1.00 


c 

0.90 

0.99 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.10 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.20 

0.14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

99 

71 

29 

44 

77 

57 

44 

74 

70 

44 

45 

29 

Got-2 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.55 

1.00 

1.00 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.45 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.50 

0.00 

0.00 

0.00 


N 

99 

54 

29 

37 

62 

47 

40 

57 

54 

64 

33 

39 
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Table IV (continuation) 


_ Station no. 

Locus - 


electromorph 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Gpi 

a 

0.00 

0.27 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.01 

0.66 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.95 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.01 

0.07 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.67 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.49 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


g 

0.00 

0.00 

0.00 

0.18 

0.19 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.82 

0.81 

0.00 

0.36 

0.32 

0.99 

0.01 

0.00 

0.00 


i 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.64 

0.01 

0.01 

0.95 

0.74 

0.99 


j 

0.00 

0.00 

0.51 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.04 

0.26 

0.01 


k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.99 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.07 

0.42 

0.47 

0.25 

0.24 

0.02 

0.40 

0.49 

0.02 

0.10 

0.44 

0.02 


N 

99 

* 

CO 

59 

91 

70 

79 

70 

119 

52 

51 

50 

45 

Hk-1 

a 

0.52 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.48 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.00 

1.00 

1.00 

1.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 


h 

0.46 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

37 

59 

38 

42 

78 

38 

28 

62 

67 

38 

67 

32 

Mdh-1 

a 

1.00 

0.90 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.00 

0.10 

1.00 

1.00 

1.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.97 

1.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 


h 

0.00 

0.16 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.05 

0.00 


N 

99 

55 

41 

61 

60 

77 

50 

54 

15 

49 

61 

39 

Mdh-2 

a 

0.00 

0.00 

0.74 

0.64 

0.48 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

1.00 

1.00 

0.26 

0.36 

0.52 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


h 

0.00 

0.00 

0.32 

0.37 

0.54 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

99 

54 

31 

40 

44 

66 

50 

58 

19 

49 

63 

51 

Me 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 


N 

99 

62 

35 

72 

85 

76 

56 

80 

22 

70 

74 

31 

Mpi 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.87 

0.13 

0.00 

0.00 

0.00 


c 

0.90 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.12 

0.87 

0.02 

0.00 

0.01 


d 

0.07 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

1.00 

1.00 

0.97 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

0.00 

0.00 

0.00 

0.31 

0.00 

0.41 


g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.67 

1.00 

0.58 


h 

0.12 

0.00 

0.00 

0.00 

0.22 

0.00 

0.00 

0.12 

0.13 

0.37 

0.00 

0.41 


N 

99 

77 

38 

53 

46 

53 

41 

64* 

35 

46 

53 

37 

Pep-2 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.92 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.48 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.08 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.52 

1.00 

1.00 


e 

1.00 

0.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


f 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.16 

0.38 

0.00 

0.00 


N 

99 

32 

27 

38 

28 

19 

48 

32 

64 

48 

18 

36 
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T Station no. 

Locus - 


electromorph 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Pep-3 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


N 

99 

32 

27 

38 

28 

19 

48 

32 

55 

49 

18 

36 

Pep-4 

a 

1.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.79 

1.00 

0.73 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.21 

0.00 

0.27 


d 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.31 

0.00 

0.33 


N 

99 

32 

CM 

38 

28 

19 

48 

32 

55 

49 

18 

30 

6-Pgd 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


N 

99 

49 

41 

36 

25 

37 

35 

48 

31 

53 

52 

39 

Pk 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.81 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.19 


h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.31 


N 

99 

34 

39 

27 

25 

21 

33 

61 

73 

38 

38 

29 


H 

0.01 

0.06 

0.04 

0.03 

0.04 

0.02 

0.04 

0.03 

0.03 

0.06 

0.03 

0.05 


numbered 1-2 (substrate 1-leucyl-glycyl-gly- 
cine), and 3-4 (substrate 1-leucyl-glycyl- 
tyrosine). The first stain (with 1-leucyl-glycyl- 
glycine as substrate) gives two zones of activity, 
the second (substrate 1-leucyl-glycyl-tyrosine) 
three. The first system in both stainings was 
identical {Pep-1). In HK and EST stains, three 
and four bands were observed, respectively. 
The Pep-1 band was not consistently scorable in 
most populations and has therefore been omit¬ 
ted. In the ADA and HK stain only the first 
(fastest) band could be scored satisfactorily; in 
the EST stain the first and the second. 

Five polymorphic loci, Apk , Hk-1 , Mpi , Pep- 
2, and Pk showed two bands in heterozygous 
individuals and were considered as monomeric 
enzymes. In six loci, Got-1 , Got-2 , Gpi , Mdh-1 , 
Mdh-2 and Pep-4 heterozygous individuals had 
three bands and these enzymes were considered 
as dimeric enzymes. No intra-population varia¬ 
tion was noted in the remaining loci (table IV). 
The Alp-2 locus was polymorphic for the two G . 
gauthieri populations from stations 6 + 7 (fig. 
2D). The subunit structure at this locus could 


not be established, but for prawns it is referred 
to as monomeric (e.g. Mulley & Latter, 1980). 

Genetic variation among populations 

Table IV shows considerable differences in 
allelic distribution among the populations 
studied. At a number of loci, the occurrence of 
specific electromorphs allows nearly complete 
discrimination between the French and the 
Iberian populations studied (viz., Apk , Est-1, 
Est-2 , Got-1, Got-2 , Pep-3 and 6-Pgd). At other 
loci, the occurrence of unique electromorphs is 
restricted to the populations from the Iberian 
peninsula: Alp-l a > c , Alp-2 a , Est-l a > b , Gdh b , Me a ’ d 
and Pkf*. G. ibericus from station 9 (tables I and 
IV; fig. ID) deviates from this general pattern 
at the Est-2 , Got-1 and Hk-1 loci. The 
monomorphic Est-2 and Got-1 loci of G. ibericus 
from station 9 are fixed for electromorphs that 
characterize the French forms. Also, in this 
sample electromorphs occur ( Ada-l a , Hk-l c and 
Got-2 a ) that have not been found in any of the 
other populations. 
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Differences between numbers of loci that are 
discriminative or diagnostic (i.e. there is no 
overlap of the allelic complements at the 0.99 
level; see Ayala & Powell, 1972) are sum¬ 
marized in table V. From the data in tables VI 
and VII, respectively, an UPGMA den¬ 
drogram and a distance Wagner network have 
been generated (figs. 3, 4, 6). From the data in 
table VI mean inter- and intra-area values of / 
can be calculated. These values are shown in 
table VIII. Samples from French and Iberian 
populations are differentiated by a high number 
of diagnostic loci (table V) and low inter-area 
values of / (table VIII). The UPGMA den¬ 
drogram (fig. 3) and the distance Wagner net¬ 
work (figs. 4, 6) show a high level of divergence 
between French and Iberian population sam¬ 
ples. Interspecific differences among French 


populations are much smaller than interspecific 
or in some cases even intraspecific differences in 
the Iberian peninsula (tables V and VIII). The 
UPGMA dendrogram and the distance Wagner 
network show clearly the higher levels of 
divergence among the Spanish populations 
compared to the French ones. As expected, 
highest /-values are found among intra-areal 
populations of G. gauthieri. Inter-areal intra¬ 
specific /-values of G. gauthieri populations are 
much lower and exhibit no overlap with intra¬ 
areal intra-specific values within this species. 
Moreover, inter-areal /-values of G. gauthieri 
populations are lower than intra-areal inter¬ 
specific values of / among populations of G. 
gauthieri and G. ibericus . There is also a con¬ 
siderable divergence between the two popula¬ 
tions of G. ibericus (stations 8 and 9: 1= 0.70). 


I- 

1.40 


—1-1-1- 

1.26 1.12 


—h- 
0.98 


—I-1-1— 

0.84 0.79 


—I-1-h- 

0.56 0.42 




€ 


1. G. stupendus 

2. G. or in os 

3- G. cf. fossa rum 

4. G.cf. fossa rum 

5. G.cf. fossa rum 

6. G. gauthieri 
G. gauthieri 
G . ibericus 
G. gauthieri 
G. gauthieri 

12. G. gauthieri 
9. G. ibericus 


—i - 1 - 1 - 1 - 1 

0.28 0.14 0.00 


Fig. 3. UPGMA dendrogram of Nei’s genetic distance (D) based upon 21 enzyme loci. 


1 -8. G. ibericus 

1 -9. G. ibericus 

I-1-1-1-f-1-1-1-1-1-1— 

0.00 0.06 0.11 0.17 0.22 0.28 


-i.o. szupenaus 

-2. G. or in os 


-a 


c 


-6. G. gauthieri 


■3. G. cf. fossarum 
4. G. cf. fossarum 
•5. G. cf. fossarum 


*7. G. gauthieri 


_J I_ 

I-17 ( 


10. <3 .gauthieri 
•11. G. gauthieri 


12. G. gauthieri 


—I-1-1-1-1-1-1-1-1 

0.33 0.39 0.44 0.50 0.55 


Fig. 4. Distance Wagner dendrogram of Rogers’ genetic distance rooted at the midpoint of the longest path (Farris, 
1972), based upon 21 enzyme loci. 
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Table V 


Number of diagnostic loci among samples of the populations studied (for further explanations see text). 


Station no. 
species 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 G. stupendus 

- 












2 G. orinos 

2 

- 











3 G. cf. fossarum 

5 

4 

- 










4 G. cf. fossarum 

5 

4 

1 

- 









5 G. cf. fossarum 

5 

4 

1 

- 

- 








6 G. gauthieri 

14 

17 

14 

15 

15 

- 







1 G. gauthieri 

14 

17 

15 

14 

14 

1 

- 






8 G. ibericus 

11 

14 

12 

11 

11 

5 

4 

- 





9 G. ibericus 

12 

14 

13 

12 

12 

9 

8 

4 

- 




10 G. gauthieri 

12 

15 

14 

15 

14 

8 

7 

12 

- 

- 



11 G. gauthieri 

13 

15 

14 

15 

14 

8 

7 

7 

11 

- 

- 


12 G. gauthieri 

12 

15 

14 

15 

14 

8 

7 

8 

12 

- 

- 

- 


Table VI 


Matrix of genetic similarity and distance coefficients. Below diagonal: Nei’s (1972) genetic distance. Above diagonal: 
Nei’s (1972) genetic identity. 


Station no. 
species 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 G. stupendus 

_ 

0.80 

0.74 

0.74 

0.75 

0.35 

0.36 

0.38 

0.39 

0.32 

0.32 

0.27 

2 G. orinos 

0.22 

- 

0.70 

0.71 

0.72 

0.27 

0.28 

0.30 

0.30 

0.25 

0.25 

0.21 

3 G. cf. fossarum 

0.31 

0.35 

- 

0.97 

0.97 

0.24 

0.25 

0.31 

0.36 

0.22 

0.22 

0.18 

4 G. cf. fossarum 

0.30 

0.35 

0.03 

- 

1.00 

0.24 

0.27 

0.33 

0.41 

0.22 

0.22 

0.18 

5 G. cf. fossarum 

0.28 

0.33 

0.03 

0.00 

- 

0.25 

0.27 

0.34 

0.41 

0.23 

0.23 

0.19 

6 G. gauthieri 

1.05 

1.30 

1.43 

1.41 

1.38 

- 

0.96 

0.76 

0.50 

0.57 

0.57 

0.60 

1 G. gauthieri 

1.03 

1.28 

1.40 

1.33 

1.30 

0.04 

- 

0.78 

0.53 

0.61 

0.60 

0.55 

8 G. ibericus 

0.97 

1.21 

1.17 

1.12 

1.09 

0.27 

0.25 

- 

0.70 

0.50 

0.59 

0.63 

9 G. ibericus 

0.95 

1.21 

1.02 

0.90 

0.88 

0.69 

0.63 

0.36 

- 

0.37 

0.37 

0.42 

10 G. gauthieri 

1.14 

1.38 

1.52 

1.51 

1.47 

0.57 

0.50 

0.53 

0.99 

- 

0.98 

0.95 

11 G. gauthieri 

1.14 

1.41 

1.52 

1.53 

1.49 

0.57 

0.51 

0.53 

0.99 

0.02 

- 

0.95 

12 G. gauthieri 

1.30 

1.57 

1.74 

1.72 

1.70 

0.51 

0.44 

0.47 

0.90 

0.05 

0.05 

- 


Main clusters in the UPGMA dendrogram and 
the distance Wagner network coincide with 
major geographical areas, although one con¬ 
spicuous difference is the placement of G. 
ibericus from station 8. According to UPGMA, 
G. ibericus from station 8 has been clustered with 
G. gauthieri from stations 6 and 7 in the same 
geographic area. According to the Wagner net¬ 
work, G. ibericus from station 9 appears more 
related to G. ibericus from station 8. However, 


the sample from station 9 is not clustered closely 
to the morphologically identical G. ibericus from 
station 8 nor with other populations from sta¬ 
tions in the same area. 

The 30 shortest (“best”) trees generated by 
the Jelly program were almost identical (e.g. 
same length, symmetry, consensus, and 
topology). The network of fig. 5 has been 
selected for being optimal. The topology of the 
Jelly network (fig. 5) is rather similar to the 
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Table VII 


Matrix of Rogers’ (1972) genetic distance. 


Station no. species 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 G. stupendus 

- 












2 G. orinos 

0.22 

- 











3 G. cf. fossarum 

0.28 

0.31 

- 










4 G. cf. fossarum 

0.28 

0.31 

0.04 

- 









5 G. cf. fossarum 

0.27 

0.30 

0.05 

0.01 

- 








6 G. gauthieri 

0.68 

0.73 

0.76 

0.76 

0.75 

- 







7 G. gauthieri 

0.65 

0.72 

0.75 

0.74 

0.73 

0.05 

- 






8 G. ibericus 

0.63 

0.70 

0.69 

0.68 

0.67 

0.25 

0.23 

- 





9 G. ibericus 

0.61 

0.69 

0.64 

0.60 

0.59 

0.51 

0.49 

0.32 

- 




10 G. gauthieri 

0.68 

0.74 

0.77 

0.77 

0.76 

0.45 

0.41 

0.42 

0.63 

- 



11 G. gauthieri 

0.68 

0.75 

0.78 

0.78 

0.77 

0.44 

0.41 

0.42 

0.63 

0.06 

- 


12 G. gauthieri 

0.73 

0.78 

0.82 

0.81 

0.81 

0.42 

0.38 

0.39 

0.60 

0.07 

0.08 

- 


Table VIII 


Mean intra- and inter-area values of the genetic identity (/) averaged by species, based upon 21 enzyme loci. G. orinos 
and G. stupendus have been put together. Inter-areal intra-specific values in bold face; range / m i n '^max * n italics. 


Station no. 

G. orinos 

G. stupendus 

1, 2 

G. cf. fossarum 

3, 4, 5 

G. gauthieri 

6, 7 

G. ibericus 

8 

G. ibericus 

9 

G. gauthieri 

10, 11, 12 

1, 2 

0.81 






3, 4, 5 

0.73 

0.98 






(0.70-0.75) 

(0.97-0.99) 





6, 7 

0.32 

0.25 

0.96 





(0.27-0.36) 

(0.24-0.27) 





8 

0.34 

0.32 

0.77 

- 




(0.30-0.38) 

(0.31-0.34) 

(0.77-0.78) 




9 

0.34 

0.39 

0.52 

0.70 

- 



(0.30-0.39) 

(0.36-0.41) 

(0.50-0.53) 

- 



10, 11, 12 

0.27 

0.21 

0.60 

0.60 

0.38 

0.96 


(0.21-0.32) 

(0.18-0.23) 

(0.56-0.65) 

(0.59-0.63) 

(0.37-0.41) 

(0.95-0.98) 



Fig. 5. Jelly “character state” network of 21 loci, based upon Rogers’ (1984) HAP algoritm (1-12 = OTU 1-12; 
13-22 = HTU 13-22). 
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2 

Fig. 6. The distance Wagner dendrogram of fig. 4 represented as an unrooted network. 

UPGMA dendrogram (fig. 3) and congruent Cross-breeding experiments 
with the distance Wagner network (figs. 4, 6). 

Electromorph frequencies of the bifurcation The results of the cross-breeding experiments 
points (hypothetical taxonomic units or are listed in table X. Experimental conditions 
HTU’s) are given in table IX. Hence, these fre- were not always optimal. Limitations on the 
quencies allow an estimation of the contribution availability of similar age groups resulted in 
of each locus to overall changes. “Character crosses between populations in different phases 
state” transformations for the loci studied are of their life cycle. Thus, in some cases, old 
indicated in table IX. specimens in the last stage of their life had to be 

Table IX 


Electromorph frequencies of hypothetical taxonomic units (HTU’s) in the Jelly network (fig. 5) at 20 enzyme loci. 
HTU’s are numbered as in fig. 5; a = fastest elecromorph; dominant electromorph(s) in bold face; * = electromorph 
p resent, frequency <0.005. 


I .nriis 

lorph 






HTU no. 





electron 

16 

14 

13 

22 

20 

15 

17 

18 

19 

21 

Ada-1 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

Alp-1 

a 

0.00 

0.00 

0.33 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

0.00 

0.00 

0.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 


c 

1.00 

1.00 

0.33 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Alp-2 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.70 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.30 

Apk 

a 

0.96 

0.97 

0.99 

0.99 

0.99 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

0.04 

0.03 

0.01 

0.01 

0.01 

1.00 

1.00 

1.00 

1.00 

0.00 

Est-1 

a 

1.00 

1.00 

0 . 00 * 

0 . 00 * 

0 . 00 * 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.99 

0.99 

0.99 

0.00 

0.00 

0.00 

0.00 

1.00 


c 

0.00 

0.00 

0.01 

0.01 

0.01 

1.00 

1.00 

1.00 

1.00 

0.00 

Est-2 

a 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

Gdh 

a 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


b 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Got-1 

a 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.02 

1.00 

1.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 


c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.90 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.08 

0.00 

0.00 

0.00 

Got-2 

a 

1.00 

1.00 

1.00 

1.00 

0.60 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

0.13 

1.00 

1.00 

1.00 

1.00 

0.00 


c 

0.00 

0.00 

0.00 

0.00 

0.27 

0.00 

0.00 

0.00 

0.00 

0.00 
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Table IX (continuation) 


HTU no. 

Locus - 


electromorph 

16 

14 

13 

22 

20 

15 

17 

18 

19 

21 

Gpi 

a 

0.00* 

0.00* 

0.01 

0.01 

0.02 

0.02 

0.23 

0.01 

0.01 

0.01 


b 

0.00* 

0.00* 

0.01 

0.07 

0.15 

0.15 

0.33 

0.05 

0.05 

0.01 


c 

0.00 

0.00 

0.00* 

0.00* 

0.00* 

0.00* 

0.02 

0.00* 

0.00* 

0.00 


d 

0.05 

0.05 

0.13 

0.26 

0.13 

0.13 

0.07 

0.04 

0.04 

0.13 


e 

0.00* 

0.00* 

0.00* 

0.01 

0.01 

0.01 

0.01 

0.00* 

0.00* 

0.00* 


f 

0.00 

0.00 

0.00* 

0.02 

0.03 

0.03 

0.02 

0.14 

0.14 

0.00* 


g 

0.05 

0.05 

0.23 

0.40 

0.55 

0.55 

0.26 

0.72 

0.72 

0.23 


h 

0.78 

0.78 

0.50 

0.18 

0.09 

0.09 

0.04 

0.03 

0.03 

0.50 


i 

0.12 

0.12 

0.01 

0.01 

0.01 

0.01 

0.01 

0.00* 

0.00* 

0.01 


j 

0.00 

0.00 

0.11 

0.04 

0.02 

0.02 

0.01 

0.01 

0.01 

0.11 


k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Hk-1 

a 

0.00 

0.00 

0.00 

0.00 

0.33 

0.62 

0.62 

1.00 

1.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

0.40 

0.28 

0.28 

0.00 

0.00 

1.00 


c 

0.00 

0.00 

0.00 

0.00 

0.27 

0.10 

0.10 

0.00 

0.00 

0.00 

Mdh-1 

a 

0.00 

0.00 

0.02 

0.05 

0.05 

0.05 

0.90 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.33 

0.87 

0.87 

0.87 

0.10 

1.00 

1.00 

0.00 


c 

0.00 

0.00 

0.33 

0.05 

0.05 

0.05 

0.00 

0.00 

0.00 

1.00 


d 

0.98 

0.98 

0.32 

0.03 

0.03 

0.03 

0.00 

0.00 

0.00 

0.00 


e 

0.02 

0.02 

0.00* 

0.00* 

0.00* 

0.00 

0.00 

0.00 

0.00 

0.00 

Mdh-2 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.48 

0.64 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.52 

0.36 

1.00 

Me 

a 

0.00 

0.00 

0.32 

0.60 

0.60 

0.00 

0.00 

0.00 

0.00 

0.00 


b 

0.00 

0.00 

0.45 

0.34 

0.34 

1.00 

1.00 

1.00 

1.00 

1.00 


c 

1.00 

1.00 

0.23 

0.06 

0.06 

0.00 

0.00 

0.00 

0.00 

0.00 

Mpi 

a 

0.00 

0.02 

0.23 

0.13 

0.03 

0.01 

0.01 

0.00* 

0.00 

1.00 


b 

0.00 

0.01 

0.16 

0.25 

0.10 

0.05 

0.05 

0.01 

0.00 

0.00 


c 

0.02 

0.03 

0.17 

0.27 

0.50 

0.48 

0.48 

0.13 

0.00 

0.00 


d 

0.00 

0.00* 

0.02 

0.04 

0.08 

0.10 

0.10 

0.03 

0.00 

0.00 


e 

0.00 

0.00* 

0.06 

0.09 

0.18 

0.26 

0.26 

0.56 

1.00 

0.00 


f 

0.30 

0.35 

0.15 

0.10 

0.08 

0.09 

0.09 

0.26 

0.00 

0.00 


g 

0.68 

0.59 

0.21 

0.12 

0.03 

0.01 

0.01 

0.00* 

0.00 

0.00 

Pep-2 

a 

0.02 

0.02 

0.10 

0.10 

0.34 

0.14 

0.14 

0.00 

0.00 

0.00 


b 

0.13 

0.13 

0.59 

0.59 

0.26 

0.06 

0.06 

0.00 

0.00 

1.00 


c 

0.06 

0.06 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.78 

0.78 

0.24 

0.24 

0.11 

0.02 

0.02 

0.00 

0.00 

0.00 


e 

0.01 

0.01 

0.06 

0.06 

0.29 

0.76 

0.76 

1.00 

1.00 

0.00 


f 

0.00* 

0.00* 

0.01 

0.01 

0.00 

0.02 

0.02 

0.00 

0.00 

0.00 

Pep-3 

a 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

0.00 

0.00 

1.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

Pep-4 

a 

0.90 

0.90 

0.99 

0.99 

0.99 

0.69 

0.69 

0.00 

0.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.09 

0.09 

0.00* 

0.00* 

0.00* 

0.00* 

0.00* 

0.00 

0.00 

0.00 


d 

0.00* 

0.00* 

0.00* 

0.00* 

0.00* 

0.08 

0.08 

0.00 

0.00 

0.00 


e 

0.01 

0.01 

0.01 

0.01 

0.01 

0.23 

0.23 

1.00 

1.00 

0.00 

6-Pgd 

a 

0.00 ‘ 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

1.00 

1.00 

Pk 

a 

0.00 

0.81 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

1.00 


b 

1.00 

0.19 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.00 
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Table X 


Results of cross-breeding experiments between populations of G. gauthieri and G. ibericus: g 7 = G. gauthieri from station 
7; gjo = gauthieri from station 10; ig= G. ibericus from station 8; i 9 = G. ibericus from station 9. 


Cross 

9 C r 

N of era 

N of 99 

N of 

ovigerous 99 

N of ovigerous 

99 with offspring 

F2 

obtained 

S7 x §7 

50 

50 

23 

5 

+ 

S10 x gio 

75 

54 

24 

7 

+ 

h X *8 

53 

45 

31 

9 

+ 

ig X ig 

100 

100 

69 

27 

+ 

glO x g7 

94 

93 

33 

11 

— 

g7 X g 10 

39 

26 

12 

7 

+ 

*'<9 X gy 

50 

50 

41 

— 

— 

g7 x ig 

76 

61 

30 

9 

+ 

ig X ig 

50 

16 

40 

8 

+ 

ig x ig 

78 

78 

69 

27 

+ 


crossed with young almost immature specimens 
from other populations. 

All crosses but one gave Fj offspring. Cross 
ig x gy (females G. ibericus from station 8 and 
males G. gauthieri from station 7) failed to yield 
offspring. Except for the reference crosses, F 2 
in very low numbers were obtained from 
crosses gyxgio (females G. gauthieri from sta¬ 
tion 7 and males G. gauthieri from station 10), 
ig x ig (females G. ibericus from station 9 and 
males G. ibericus from station 8), ig x ig (females 
G. ibericus from station 8 and males G. ibericus 
from station 9) and gy x ig (females G. gauthieri 
from station 7 and males G. ibericus from station 
8 ). 

DISCUSSION AND CONCLUSIONS 
Genetic diversity 

Mean /-values describing various taxonomic 
levels are given by several authors (Ayala et al., 
1973; Avise & Smith, 1977; Ferguson, 1980; 
Thorpe, 1982; Menken & Ulenberg, 1987; see 
table XI). Inter-areal values of 7 between 
French and Iberian populations separated by 
the Pyrenean Mountain chain are remarkably 
low (7=0.18-0.41, mean: 0.28; 7) = 0.88-1.75, 
mean: 1.30). Referring to table XI, these 
values locate at the generic level. Similar results 


were reported by Siegismund et al. (1985) and 
Skadsheim & Siegismund (1986) who stated that 
“within the genus Gammarus , genetic identities 
in species pairs are as low as those between 
Gammarus spp. and species in other genera”. 

Intra-specific values of 7 among G. gauthieri 
populations from the Valladolid area and 
Cuenca area (7=0.56-0.65; D = 0.57-0.44), 
and G. ibericus populations from stations 8 and 
9 (7=0.70; D = 0.36) are low for conspecific 
populations. The Valladolid area and the Ser- 
rania de Cuenca are separated by the Sierra de 
Gredos and the Sierra de Guadarrama. Stations 
8 and 9, although geographically much closer, 
belong (see table I) to the Atlantic Ocean and 
Mediterranean drainage systems, respectively, 
separated by the Serrania de Cuenca. These 
mountain chains have also been recognized by 
Lop (1987) as barriers preventing genetic 
exchange. According to table XI, these values 
are below the expected subspecies level and 
range about what is commonly attributed to 
sibling species, and, to a lesser extent, con¬ 
generic species. 

Intra-specific values of D in G. gauthieri range 
from 0.02-0.57 (7=0.56-0.98). Lop (1987), in 
his study of the Echinogammarus berilloni group 
from the Iberian Peninsula, found intra-specific 
values of D ranging from 0.05-0.26 (7=0.77- 
0.95). Although these values diverge less than 
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Table XI 

Values of genetic identity on various taxonomic levels based upon literature data: 1 - Ayala et al. (1973); 2 - Avise & 
Smith (1977); 3 = Thorpe (1982); 4 = Ferguson (1980); 5 « Menken & Ulenberg (1987); /= genetic identity (Nei, 1972); 
S.E. = standard error. 


Taxonomic level 

1 

/ S.E. 

2 

/ S.E. 

3 

I S.E. 

4 

/ 

5 

Range of I 

Populations 

0.970 ±0.006 

0.980 ±0.004 

0.960 ±0.032 

— 

0.78-1.00 

Subspecies 

0.795 ±0.013 

0.843 ±0.013 

— 

— 

0.80-0.98 

Semispecies 

0.873 ±0.012 

— 

— 

— 

— 

Sibling species 

0.517 ±0.024 

— 

— 

0.626 

0.00-0.98 

Species 

0.352 ±0.0023 

0.544 ±0.020 

0.540 ±0.168 

0.567 

0.00-0.90 

Genera 

— 

0.294 ±0.020 

0.273 ±0.107 

0.362 

0.00-0.60 


in G. gauthieri , it can be concluded that gam- 
mar ids from the Iberian Peninsula tend to be 
genetically much more differentiated at the 
species level than species groups reported in the 
literature, or even Gammarus spp. in France: the 
/-value for G . orinos and G. stupendus is 0.80 
(.D = 0.22) and other data suggest that lower 
values are not likely to be found (Scheepmaker, 
in prep.). There is a discrepancy between mor¬ 
phological and allozymic differentiation. Evolu¬ 
tion of morphological characters may not be 
constant. If particular morphological types tend 
to be conserved, the relation between mor¬ 
phological and genetic divergence will decrease. 

Morphologically, G . cf. fossarum is related to 
G. orinos rather than to G. fossarum Koch. From 
this study it may appear that G. cf. fossarum 
should be considered as a separate species. Fur¬ 
ther experiments showed that G. cf. fossarum 
from the Pyrenees and G. fossarum are not con- 
specific (Scheepmaker, in prep.). In these 
experiments however, some populations of G. 
cf. fossarum and G . orinos showed a considerable 
overlap in genetic variation. Thus, the tax¬ 
onomic status of G. cf. fossarum remains provi¬ 
sionally unsolved and will be treated in a future 
study. 

Genetic divergence and evolutionary rates 

Several authors have attempted to explain dif¬ 
ferences in genetic variation of enzymes within 
and among populations and species (e.g. 


metabolic function: Johnson, 1976; Selander, 
1976; for Crustacea, Hedgecock et al., 1982; 
molecular structure and subunit size: Koehn & 
Eanes, 1978). Evidence for certain loci being 
more conservative than others has been dis¬ 
cussed by Avise (1975), Sarich (1977), Thorpe 
(1982) and Menken & Ulenberg (1987). Sarich 
(1977) distinguished “fast” and “slow” evolv¬ 
ing loci. According to Menken & Ulenberg 
(1987), the ability to differentiate among the 
investigated species may depend on different 
evolutionary rates. Striking is, in this study, 
that a number of loci appear to be fixed for 
specific electromorphs in many samples in 
which the alternative allelic complement is rare, 
or absent (this was also found in other studies 
on amphipods, e.g. Siegismund et al., 1985; 
Bulnheim, 1985; Bulnheim & Scholl, 1981, 
1986; Scheepmaker, unpubl.). This may be 
consequent to a local polymorphism fixed for 
the mutant electromorph by drift or directional 
selection. In either case, it may be justified to 
calibrate the molecular clock for a given set of 
loci in a group of closely related species, as long 
as the evolutionary rates for the individual loci 
remain constant. 

These “fixed” loci are invariable or almost 
invariable across the populations of a species or 
species group over long distances. For instance, 
the Lap locus is invariable for all members 
studied of the G. pulex group (sensu Karaman & 
Pinkster, 1977) ranging from Holland to cen¬ 
tral Spain. Within the same group, the 6-Pgd 
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and Pep-3 loci are fixed for one electromorph in 
all members from the Iberian Peninsula (and 
thus diagnostic with respect to those north of 
the Pyrenees). All the other members, from the 
Pyrenees north up to Holland, are fixed 
invariably for the same alternative elec¬ 
tromorph. At the species level, electromorphs 
are sometimes invariable over long distances 
across most members, but diagnostic for one or 
two of them (e.g. G. fossarum, G. orinos, G. cf. 
fossarum and G. stupendus versus G. wautieri and 
G. pulex at the Pk locus; Scheepmaker, in 
prep.). Occasionally, this was also observed at 
the intra-specific level (viz. Ada-1 and Hk-l c at 
station 9; Pk at station 12). 

In the present study, such loci seem to be rep¬ 
resented by Apk, Ada-1, Alp-1, Alp-2, Est-1, Est - 
2, Gdh, Hk-1, Lap, Pep-3, 6-Pgd, and Pk. Alter¬ 
natively, at both inter- en intra-specific level, 
loci such as Gpi, Mpi, Pep-2, Pep-4, Got-1, Got-2 , 
Mdh-1, and Mdh-2 may be highly variable. The 
remaining loci may range somewhere in 
between. This bimodality in levels of polymor¬ 
phism (and possibly in rates of evolutionary 
change) has been reported by several authors 
(reviewed by Thorpe, 1982 and Menken & 
Ulenberg, 1987). 

When calculating distance matrices (and the 
resulting dendrograms) no distinction is made 
between slow and fast evolving loci. However, 
since the Wagner and UPGMA dendrograms 
and the Jelly network are similar in overall view 
(figs. 3, 4, 6), major splits and path-lengths can 
be brought back to ‘ ‘character state” changes 
in the Jelly network (fig. 4, table IX). Employ¬ 
ing the Jelly computer program, an estimate of 
the amount of change and the proportion of fast 
and slow loci involved can be made. 

Major splits in the Jelly network are 
generated by loci that either exhibit a low level 
of polymorphism, or are monomorphic (e.g. 
Apk, Alp-1, Alp-2, Est-1, Est-2 , Gdh, Hk-1 , 
6-Pgd, Pk; see table IX). The contribution of 
highly variable loci such as Gpi, Mpi, Pep-2, and 
Pep-4 is reduced. Alternatively, these loci 
(reflecting local variation in conspecific popula¬ 
tions) determine largely the topology of the net¬ 
work at the OTU level (e.g. G. cf. fossarum from 


stations 3, 4, 5; G. gauthieri from stations 6 and 
7, and stations 10, 11, 12). 

Cross-breeding experiments 

In cross-breeding experiments, French and 
Iberian species have proved to be completely 
infertile among each other (e.g. G. ibericus vs. 
G. orinos, Pinkster & Scholl, 1984; G. gauthieri 
from Spain vs. G. wautieri and G. fossarum, 
Goedmakers & Roux, 1975). However, the 
results of the cross-breeding experiments in this 
study (table X) do not contribute to the resolu¬ 
tion of the phylogenetic relationships or the tax¬ 
onomic status of genetically well differentiated, 
but morphologically similar allopatric forms of 
G. gauthieri and G. ibericus from the Iberian 
Peninsula. 

Despite of the considerable amount of genetic 
divergence between allopatric populations of 
both G. gauthieri and G. ibericus, complete 
reproductive isolation has not yet been 
achieved. However, it can also be concluded 
that these conspecific allopatric populations are 
not completely interfertile. Another remarkable 
result is provided by cross gy x ig (females G. 
gauthieri from station 7 and males G. ibericus 
from station 8), yielding F2 offspring. These 
populations are both genetically and mor¬ 
phologically clearly distinct. 

It could be argued that morphologically iden¬ 
tical and only partially reproductively isolated 
forms should be considered conspecific. How¬ 
ever, in the case of genetically and mor¬ 
phologically well-differentiated forms such as 
G. gauthieri from station 7 and G. ibericus from 
station 8, which still prove to be at least par¬ 
tially interfertile, reproductive isolation as 
criterion in taxonomy becomes heavily 
undermined. 

Furthermore, it is interesting that both mor¬ 
phologically and genetically much less differen¬ 
tiated species in France are often reproductively 
isolated completely (e.g. G. pulex v s. G. fossarum, 
Wautier & Roux, 1959; Roux, 1971a; G. 
wautieri vs. G. pulex and G. fossarum, Roux, 
1971a; G. stupendus vs. G. orinos and G. fossarum, 
Pinkster, 1983; Scheepmaker, 1987; G. orinos 
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vs. G. cf. fossarum and G. Joss arum, Scheep- 
maker, in prep.). For the present study, cross¬ 
breeding experiments were also performed with 
G. cf. fossarum from station 3 (fig. 1B) and G. 
ibericus from station 9 showing electromorphs 
that characterize the French forms. The 
technical difficulties mentioned above, and 
increased mortality in one of the populations 
involved, caused these experiments to be 
unsuccessful. 

Roux (1971b) demonstrated that sympatric 
occurrence enhances reproductive isolation. 
This suggests that reproductive isolation may 
not be related to the degree of genetic differen¬ 
tiation, but rather to the mode of speciation 
(e.g. allopatric vs. sympatric). 

Genetic distance and divergence time 

Nei’s D claims to be a measure for genetic 
divergence time (Nei, 1975). However, it is 
obvious that the divergence time thus obtained 
is strongly dependent on the number and the 
proportions of fast and slow loci involved and 
the organisms studied. Moreover, the same loci 
appear to contribute to a different extent to D - 
values in comparable stages of evolutionary 
divergence depending on the taxonomic group 
studied (De Matthaeis et al., 1983). These 
reasons might explain in part the discordance of 
approximations of D that have been made in 
the past: proposed values of one unit D range 
from 5 to 20 mY (Thorpe, 1982; Sbordoni et 
al., 1980; Bert, 1986). 

Sbordoni (1982), who reviewed speciation in 
cave animals, concluded that troglobytes, being 
confined to caves and not having the oppor¬ 
tunity to disperse over long distances, are 
excellent paleogeographic monitors and their 
present distribution generally reflects particular 
paleogeographic events which have been dated 
by geologists. Several authors (Sbordoni et al., 
1980; Caccone et al., 1982; De Matthaeis et al., 
1982; Cobolli Sbordoni et al., 1987; Caccone et 
al., 1986) studied divergence time in cave¬ 
dwelling organisms in the Mediterranean 
region. These authors reported a close agree¬ 
ment of paleogeographic and paleoclima- 


tological evidence and divergence time ac¬ 
cording to Nei (1975) (viz. a one unit D- 
value of 5 mY). Thus, from these studies, with 
in most cases related groups of organisms, it 
appears that Nei’s (1975) estimate of 
divergence time provides a reasonable approx¬ 
imation. Therefore, to check if there is any 
agreement between ranges of divergence time 
and paleogeographical and paleoclimatological 
evidence, we calculated divergence time 
according to Nei (1975). 

Two alternative hypotheses can be proposed 
to account for the present-day divergence of 
French and Iberian members of the G. pulex 
group. According to the first, French popula¬ 
tions of G. orinos , G . stupendus and G. ci. fossarum 
and Iberian populations of G. gauthieri and G. 
ibericus may have originated from the same 
stock and separated in the past by some barrier. 
Evidence for the Pyrenean mountain chain as a 
barrier preventing gene flow is provided by 
high D- values between species from France and 
the Iberian Peninsula. However, calculated 
divergence time should be in the same order as 
geological data on the uplifting of the Pyrenees. 
The Alpine orogenesis in the Pyrenean belt 
started in the Late Cretaceous. About that 
period the Iberian Peninsula was intermittently 
separated from the European continent by 
marine transgressions, and during the Santo- 
nian, the Alpine sea was connected with the 
Atlantic (Plaziat, 1981). Although folding first 
occurred at the end of the Cretaceous, major 
folding occurred at the end of the Oligocene 
and vertical movements resulting into the 
present-day relief may be largely of Pleistocene 
age, thus less than 2 mY ago (Mattauer & 
Henry, 1974). It is not clear in which period the 
Pyrenees became a barrier to gene flow. 
Minimum and maximum Z)-values according 
to Nei (1975) among the Spanish and French 
populations studied range from 8.7 to 4.4 mY 
(Upper Miocene-Middle Pliocene; mean 6.5 
mY). 

If the Pyrenees formed a barrier to gene flow, 
the fall in sea level during the Messinian 
salinity crisis (Hsu, 1972; Hsu et al., 1977; 
Rogl & Steininger, 1983) possibly provided a 
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temporary migration pathway along the emerg¬ 
ing coasted plains bordering the eastern bound¬ 
aries of southern France and Spain; this may 
also have happened during the Pleistocene 
glaciations and the subsequent falls in sea level. 
Supporting evidence may be provided by the 
‘‘French’* electromorphs Esi-2 b , Got-l a and 
Got-F in G. ibericus from station 9. Moreover, 
the presence of this electrophoretically distinct 
form coincides with the southernmost distribu¬ 
tion limit in this area of the western European 
Echinogammarus berilloni group (Lop, 1987). 
During the Messinian, G. gauthieri was probably 
able to disperse to North Africa. 

Although divergence data and geological 
evidence do not conflict, they are hardly infor¬ 
mative and leave the way open to a second 
hypothesis. According to this hypothesis, both 
French and Iberian populations originated 
from the same stock, probably in the Balkan 
and Asia Minor, areas considered as “centre of 
origin” of the G. pulex group because of its 
species richness (Karaman & Pinkster, 1987). 
According to this theory, (ancestors of) West 
European species (viz. the French species) 
moved along the northern side of the Mediter¬ 
ranean to the west, whereas the Iberian species 
originated from ancestors that migrated along 
the North African side of the Mediterranean 
after the junction of the Eurasian and African 
plates in the Middle Miocene (about 10 mY 
ago). During the Messinian salinity crisis (6-5 
mY ago) these forms where able to migrate 
from North Africa to Spain. Supporting 
evidence for this theory is provided by the 
distribution areas of G. gauthieri and the 
Echinogammarus simoni group, which both have 
their largest extension in North Africa 
(Karaman & Pinkster, 1977). 

The hypothesis proposing that the ancestors 
of G. gauthieri and G. ibericus originated from 
Europe and not from Africa also agrees with the 
“paleo-ecological principle” proposed by 
Golikov & Tzvetkova (1972). Because of the 
improbability of a major adaptive shift of a 
species group as a whole to different ecological 
conditions, this principle presumes that the 
major part of the contemporary distribution 


area of a species group is still limited primarily 
to ecological conditions that are similar to those 
in which it originated. For gammarids, these 
particular conditions are found in the more 
temperate regions. During the Middle 
Miocene, perhaps except at higher altitudes, 
temperate conditions were not yet prevailing in 
North Africa. Hence, according to this 
hypothesis, the extension of gammarids is most 
likely to represent an adaptive shift of some 
species to the more extreme conditions in a 
remote area rather than evidence for a 
presumed alternative centre of origin of the 
species group. 

As a third alternative, a combination of both 
dispersion routes may account in part for the 
present-day distribution. 

Paleoclimatological aspects of inter-areal genetic 
divergence 

A general trend towards a cooler and dryer 
climate occurred during the Late Miocene 
(Maldonado, 1985; Hsii et al., 1977), while 
during the Pliocene a sequence of cool and cold 
phases started (Van der Hammen et al., 1971; 
Lop, 1987). This sequence preceded the 
Pleistocene glaciations. During the cold phases, 
major shifts in vegetation zones occurred. 
North of the Pyrenees boreal conditions may 
have prevailed whereas the Iberian Peninsula 
was characterized by a steppe and desert 
climate (Van der Hammen et al., 1971). In the 
Iberian Peninsula the semi-permanent desicca¬ 
tion of riverbeds and watercourses as a conse¬ 
quence of the dry climate may have resulted in 
long-term isolation of populations and conse¬ 
quently in increased levels of genetic 
divergence. This has also been suggested by 
Lop (1987) who found high levels of intra- 
specific divergence in Echinogammarus spp. from 
the Iberian Peninsula. 

Alternatively, under the more humid boreal 
conditions north of the Pyrenees, the oppor¬ 
tunity to disperse over long distances was main¬ 
tained. This hypothesis appears to be supported 
by Z)-values among the Spanish and French 
populations studied. The mean Z)-value among 
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Spanish populations is 0.47 with a maximum of 
1.00, which corresponds to about 2.3 and 5 
mY, respectively, and suggests a divergence of 
Pliocene age. The mean D -value of French 
populations is 0.22 (maximum 0.35), corre¬ 
sponding to 1.1 and 1.7 mY, respectively, sug¬ 
gesting a divergence of Pleistocene age. In this 
period, cold phases were likely to be responsible 
for an increased isolation among populations 
north of the Pyrenees. 

It may be concluded that the orders of 
divergence resolved and data on paleogeo- 
graphical and paleoclimatological evidence do 
not conflict. These data provide an explanation 
for differences in inter-areal genetic divergence 
of French and Iberian populations. However, 
they do neither contribute to the solution of 
dispersion routes nor to the appreciation of dif¬ 
ferences in the evolutionary rates of the loci 
involved. Electrophoretic screening of North 
African populations of G . gauthieri seems 
indispensable to resolve the former problem, 
and to potentially provide an initial calibration 
of the molecular clock for the G. pulex group for 
the given set of loci. 
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